Radio Polarization of BL Lacertae objects 
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Abstract 

In this paper, using the database of the university of Michigan Radio Astronomy Observatory (UMRAO) 
at three (4.8 GHz, 8 GHZ, and 14.5 GHz) radio frequencies, we studied the polarization properties for 47 
BL Lacertae objccts(38 radio selected BL Lacertae objects, 7 X-ray selected BL Lacertae, and two inter- 
middle objects (Mkn 421 and Mkn 501), and found that (1) The polarizations at higher radio frequency is 
higher than those at lower frequency, (2) The variability of polarization at higher radio frequency is higher 
than those at lower frequency, (3) The polarization is correlated with the radio spectral index, and (4) The 
polarization is correlated with core-dominance parameter for those objects with known core-dominance 
parameters suggesting that the relativistic beaming could explain the polarization characteristic of BL 
Lacs. 
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1. Introduction 

BL Lacertae objects are generally described as a sub- 
class of active galactic nuclei (AGNs), showing rapid and 
large amplitude variation, variable and high polarization, 
and core-dominated non-thermal continuum. Some BL 
Lacertae objects show supcrluminal motion and high en- 
ergy gamma-ray emissions ( Andruchow et al. 2005; Angel 
& Stockman 1980; Fan 2005; Fan et al. 1996; Gabuzda 
2003; Gu et al. 2006; Gupta et al. 2004; Qian et al. 2003, 
2004; Romero et al. 1995, 2002; Stickel et al. 1993; Tao 
& Qian, 2004; Wills et al. 1992; Xie et al. 2004). 

According to the surveys, BL Lacertae objects can be 
divided into radio selected BL Lacertae objects (RBLs) 
and the X-ray selected BL Lacertae objects (XBLs). From 
the spectral energy distribution (SED), one can see that 
some BL Lacertae objects show lower-frequency peaked 
SED and are called as LBLs, while some others show 
higher-frequency peaked SED and are called as HBLs (see 
Urry & Padovani 1995 for a review). Generally, XBLs cor- 
respond to HBLs while RBLs to LBL. For Mkn 421 and 
Mkn 501, they are classified as RBLs by someone and as 
XBLs by some other authors. However, from SED, they 
are both classified as HBLs. 

The observational properties of RBLs are systemati- 
cally different from those of XBLs. The latter have flatter 
spectral energy distribution from the radio through X-ray, 
a higher observed peak of the emitted power from radio 
through X-ray spectral energy distribution. Furthermore, 
RBL and XBL show different observation properties(see 
Angel & stockman 1980; Efimov, et al. 1988a,b, 2002; 
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Fan et al. 1997; Giommi et al. 1995; Junnuzi et al. 1994; 
Sambruna et al. 1996). 

In this paper, we will investigate the radio properties of 
polarization in BL Lacertae objects based on the MURAO 
data base. It is arranged as follows. In section 2, we will 
calculate the averaged polarization and the variation of 
polarization for RBLs and XBLs. In section 3, we give 
some discussion and a conclusion. 

2. Calculation and Results 

From the MURAO data base, we can get 47 BL Lacertae 
objects, 38 of them are RBLs, 7 of them are XBLs, two of 
them (Mkn 421 and Mkn 501) arc intermiddlc type. For 
those sources, we did following calculations 

2.1. Radio Spectral Index 

For each source, firstly we got its averaged flux densities 
at three radio frequencies (4.8GHz, 8.0GHz and 14.5 GHz) 
respectively, then adopted the linear regression analysis to 
the averaged flux at the three frequencies, we got the fit- 
ting result as logS*,, = — alog^ + c. Therefore, we finally 
got the spectral index, a (S v oc v~~ a ). The resulting spec- 
tral indexes are listed in Table 1, from which, we found 
that the radio spectral indexes for the 7 XBLs are in the 
range of a = -0.188 to a = 0.625 with an averaged value 
of (a) = 0.235 ±0.255 while those for the 38 RBLs are in 
the range of a = -0.601 to a = 0.757 with an averaged 
value of (a) = 0.044 ± 0.264. 

2.2. Radio Polarization 

For the RBLs and XBLs, we calculated respectively the 
averaged maximum polarization and the averaged polar- 
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ization based on the observed polarization of each source 
The detail procedure is as follows. At each frequenc; 
for each source, we chose the observed maximum polar 
ization value as the maximum polarization at the corrc 
sponding frequency, and calculated the averaged value c 
the whole observation polarization as the averaged po 
larization. Therefore, we have three maximum polariza 
tions and three averaged polarizations (at 4.8, 8.0 ani 
14.5 GHz) for each source. Then for each subclass (XBL 
and RBLs), we can calculate the averaged maximum po 
larization at frequency v (v = 4.8, 8.0 and 14.5 GHz) 
pMax = Y^pM ax /N, here P^ ax is the maximum observe, 
polarization at frequency v for the ith source in the sub 
class, and N is the source number of the subclass. N = 
and 38 for XBLs and RBLs respectively. The averaged av 
eraged polarization is also calculated as P^ ve = T,P/^ e /N 
here P^ e is the averaged observed polarization for the it] 
source in the subclass. So, we can get following statistica 
results. 

«80GHz(%))=36±29, 
( J P8 M 00GHz(%))=65±23, 

(^.Tgh z (%))=74±25, 



and 



for 7 XBLs, 



<^80GH Z (%))=7±5, 
<^8 A 00GHz(%)Hl6±7, 
(A A 4 V 5GHz(%)) = 17±7 



< J P4 M 80GHz(%)) = 17±13, 

<^8 M ooghz(%))=26±21, 

«.5GHz(%))=29±21, 



and 



(' P 4 A 80GHz(%)) 



(^OOGHzTO) 

/ pAve 
\- r 14.5GHz 



= 4±2, 
= 5 ±3, 
(%)>=6±3 



for 38 RBLs. 
2.3. Results 



For polarization, the averaged results suggest that botl 
the maximum and averaged polarizations increase wit! 
the radio frequency in the radio range of from 4.8GHz ti 
14.5GHz. In addition, within that radio frequency range 
both the averaged maximum and the averaged average* 
polarization in XBLs are higher than those in RBLs a 
the same radio frequency. 

For the polarization and the spectral index, there is i 
tendency for the radio polarization to increase with tb 
spectral index for BL Lacertae objects as shown in Fig. 1 

The core-dominance parameter is an important param 
eter calculated from the radio observation. It is defined a a 
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Fig. 1. The relation between the averaged polarization at 
8 GHz and the radio spectral index for radio selected BL 
Lacertae objects-RBL (triangles) and X-ray selected BL 
Lacertae objects— XBL (open circles) 
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Fig. 2. The relation between the averaged polarization at 
8 GHz and the core-dominance parameter for radio selected 
BL Lacertae objects-RBL (triangles) and X-ray selected BL 
Lacertae objects-XBL (open circles). Curve 1 stands for r\ = 
0.25, curve 2 stands for 77 = 0.18, curve 3 stands for rj = 0.11, 
and curve 4 stands for r) = 0.05. 

the ratio of the core flux density to the extended flux den- 
sity (Orr &; Browne 1982). For polarization and the core- 
dominance parameter, we plotted the averaged 8GHz po- 
larization against the core dominance parameter as shown 
in Fig. 2. 

3. Discussion 

It is believed that the particular properties observed 
from BL Lacertae objects are associated with the beam- 
ing effect, and the beaming model was adopted to explain 
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Fig. 3. The relation between the averaged polarization at 
8 GHz and the polarization variability-erp for radio selected 
BL Lacertae objects-RBL (triangles) and X-ray selected BL 
Laccrtae objects— XBL (open circles). The line stands for the 
best fitting result. 

both the particular observational properties and some ob- 
servational differences between RBLs and XBLs (see Fan 
et al. 1997; Fan & Xie 1996; Georganopoulos & Marscher 
1999) although the viewing angle alone can not explain 
all the difference between RBLs and XBLs (Sambruna et 
al. 1996). 

High and variable polarization is one of the character- 
istics of BL Lacertae objects. For the radio polarization, 
although the averaged value indicates that the polariza- 
tion in XBLs is higher than that in RBLs, we can get the 
conclusion that the radio polarization in XBLs is not al- 
ways higher than that in RBLs. The reasons are that 1) 
There are only 7 XBLs in our consideration, and 2) for 
the sources with the averaged polarization being higher 
than 10%, we have two RBLs but 6 XBLs. However, for 
the sources with maximum polarization being higher than 
35%, we have 9 RBLs but 3 XBLs. We think that whether 
or not the radio polarization in XBLs is higher than that 
in RBLs, more radio polarization observations are signifi- 
cant and encouraged. 

When we take into account of the polarization variation, 
we can use the standard deviation, erp as the variation 
indication. In this sense, we found that the variation is 
closely linearly correlated with the averaged polarization 
shown in Fig. 3, the line stands for the best fitting result, 
P avc (%) = (1.02 ±0.04)o- P + (1.54 ±0.32). Here we used 
the data at 8GHz. 

When we considered the relationship of polarization on 
the spectral index, we found that the radio polarization 
increases with the spectral index, a (S^ oc v~ a ). The data 
can gave a linear relation of P ~ 9.15a + 7.04. From the 
synchrotron emission, one can get that the polarization 
is proportional to a ^_~t^ 3 , which suggests that the polar- 
ization is correlated with the spectral index and that the 



polarization increases with the spectral index. Although 
our statistical result show a trend that the polarization 
increase with the radio spectral index as shown in Fig. 1 , 
the trend is largely different from the prediction by syn- 
chrotron emission. 

Polarization is found to be associated with the core- 
dominance parameter ( see Wills et al. 1992 and ref- 
erence therein) with high polarization corresponding to 
large logi?. For the present sample, one can get a plot 
of the radio polarization against the core-dominance pa- 
rameter as shown in Fig. 2. In Fig. 2, we only showed 
the averaged radio polarization at 8.0GHz and the core 
dominance parameters that are available in the literature. 

Based on the two-component beaming model Urry & 
Shafer (1984) (see also Urry & Padovani, 1995), the emis- 
sions of an AGN are from two components, namely the 
beamed and the unbeamed ones. Assuming the intrinsic 



to be a fraction / of the unbeamed 

ob q i cob 



flux of the jet, 5] 

flux, S unb , i.e., Sp = fS u „b, one can get 5°" = 5 Mrlfa + 5 J 
(1 + fS p )S un b- If one assumes that the emissions in the co- 
moving jet are also composed of two components, namely 
the polarized and the unpolarized, with the two compo- 
nents being proportional to each other, S™ = S? + Sj p , 
= vSj P - Then a relation between the observed polar- 
ization and the Doppler factor can be obtained (Fan et al. 
1997; see also Fan et al. 2001). 
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is the intrinsic polarization expressed in the 
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(3) 



1 + fSP 1 + v 

S is the Doppler factor, r\ is the ratio of the polarized to the 
unpolarized luminosity in the jets. The value of p depends 
on the shape of the emitted spectrum and the detailed 
physics of the jet (Lind & Blandford 1985), p = 3 + a is 
for a moving sphere and p = 2 + a is for the case of a 
continuous jet, a is the spectral index. 

If we use the expression for the core-dominance param- 
eter, R = f((r(l-j3cos6))-P + (r(l+l3cos6))- p )=fSP[l + 
f( \+f°ose ) P ^ thcn thc polarization can be written in thc 
form 



pob _ 



R-GU,0,P) V 
l-G(f,9,j3) + Rl + r) 

,P) = f{\- 0cos6 



(4) 



Since the observed polar- 



here, G{f,L, yj - jK^^gj 
ization is positive, then we can get the value of G(f,8,f3) 
to be smaller than the value of R. Therefore, thc value of 
G(f,8,/3) is actually very small since thc logi? is -0.91 for 
1100+772. So, we have 



P 



ob 



R 



l + Rl + ri' 
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For a given 77, a theoretical relation of polarization depend- 
ing on the core-dominance parameter can be obtained. In 
Fig. 2, we showed the polarization vs core-dominance pa- 
rameter plot. The core-dominance parameter is taken as 
the jet viewing indication, it is actually the indication of 
a beaming effect. For a strongly boosted source, the core- 
dominance parameter can be simply expressed as R = fS p 
(Ghisellini et al. 1993). In this sense, one can see clearly 
that 1) the polarization is consistent with the explanation 
that the polarization is associated with the beaming effect 
and 2) for most sources, their r\ is less than 0.25. 

In this paper, based on the MURAO data base, we con- 
sidered the radio polarization for BL Lacertae objects and 
found that the averaged polarization depends on the fre- 
quency, the spectral index and the core-dominance param- 
eter. The polarization variation is closely correlated with 
the polarization. The rclativistic beaming could explain 
the polarization characteristic of BL Lacs. 
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Table 1. A sample of 47 BL Lacertae objects 



Source 


Class 


a 


Tjave. 
r SGHz 


a p 


Data span 


logi? 


Ref 


0003-066 


RBL 


0.088 


3.27 


1.54 


1978- 


1998 






0048-097 


RBL 


-0.242 


5.4 


3.05 


1970- 


1999 


1.3 


G93 


0109+224 


XBL 


-0.188 


7.52 


5.73 


1978- 


1999 


1.62 


W92 


0215+015 


RBL 


-0.337 


3.16 


2.31 


1979- 


1999 


0.9 


W92 


0219+428 


RBL 


0.754 


2.96 


2.3 


1974- 


1999 


0.3 


G93 


0235+164 


RBL 


-0.192 


2.74 


2.13 


1974- 


1999 


2.2 


G93 


0300+470 


RBL 


-0.058 


3.38 


1.99 


1975- 


1999 


1.9 


G93 


0323+022 


XBL 


0.252 


25.02 


22.55 


1985- 


1999 






0422+004 


RBL 


-0.191 


5.85 


5.35 


1978- 


1999 


3.3 


W92 


0716+714 


RBL 


-0.271 


9.19 


8.97 


1981- 


1999 


-0.52 


W92 


0735+178 


RBL 


0.062 


2.48 


1.49 


1977- 


1999 


3.4 


G93 


0754+100 


RBL 


-0.09 


5.3 


2.62 


1978- 


1999 


1.14 


W92 


0808+019 


RBL 


-0.133 


5.98 


3.06 


1979- 


1999 


1.55 


W92 


0814+425 


RBL 


-0.035 


4.06 


2.18 


1977- 


1999 


0.9 


W92 


0818-128 


RBL 


-0.054 


4.47 


4.08 


1979- 


1999 


0.23 


W92 


0829+046 


RBL 


-0.239 


3.42 


1.78 


1978- 


1999 


0.71 


F03 


0851+202 


RBL 


-0.385 


4.78 


2.83 


1971- 


1999 


3.5 


G93 


0912+297 


XBL 


0.318 


12.06 


8.65 


1980- 


1999 


0.3 


W92 


0954+658 


RBL 


-0.022 


14.5 


16.11 


1974- 


1999 


0.53 


W92 


0957+227 


RBL 


0.757 


12.74 


12.36 


1979- 


1999 






1100+772 


XBL 


0.625 


16.17 


11.95 


1983- 


1999 


-0.91 


W92 


1101+384 


Mid-BL 


0.182 


4.94 


3.85 


1978- 


1999 


0.76 


F03 


1133+704 


XBL 


0.184 


26.41 


18.67 


1980- 


1999 






1147+245 


RBL 


0.023 


3.85 


2.72 


1979- 


1999 


1.42 


W92 


1215+303 


RBL 


-0.008 


6.22 


6.05 


1979- 


1999 


0.27 


W92 


1219+285 


RBL 


-0.065 


4.75 


3.02 


1978- 


1999 


-0.72 


F03 


1307+121 


RBL 


-0.002 


5.88 


3.28 


1978- 


1999 






1308+326 


RBL 


-0.205 


2.72 


1.56 


1976- 


1999 


1.7 


G93 


1400+162 


RBL 


0.481 


8.99 


8.51 


1978- 


1999 


-0.11 


F03 


1413+135 


RBL 


-0.601 


2.2 


2.12 


1978- 


1999 


0.55 


F03 


1418+546 


RBL 


-0.117 


3.28 


3.35 


1978- 


1999 


1.77 


W92 


1514+197 


RBL 


-0.065 


6.75 


4.57 


1978- 


1999 


2.18 


W92 


1538+149 


RBL 


-0.147 


5.95 


2.21 


1977- 


1999 


0.95 


W92 


1652+398 


Mid-BL 


0.208 


3.5 


2.03 


1977- 


1999 


0.07 


F03 


1717+178 


RBL 


-0.127 


7.26 


5.07 


1977- 


1999 


1.78 


W92 


1727+502 


XBL 


0.41 


20.9 


20.84 


1979- 


1999 


1.01 


W92 


1749+096 


RBL 


-0.473 


3.18 


2.15 


1978- 


1999 


0.5 


F03 


1749+701 


RBL 


0.051 


8.45 


6.15 


1980- 


1999 


-0.46 


F03 


1803+784 


RBL 


-0.102 


4.27 


2.23 


1981- 


1999 


0.42 


F03 


1807+698 


RBL 


0.048 


4.54 


1.66 


1979- 


1999 


0.68 


F03 


1823+568 


RBL 


-0.099 


6.77 


2.52 


1981- 


1999 






2007+777 


RBL 


-0.091 


5.22 


2.5 


1981- 


1999 


1.36 


G93 


2032+107 


RBL 


0.116 


4.55 


4.6 


1978- 


1999 


1.61 


W92 


2131-021 


RBL 


0.065 


2.34 


1.66 


1974- 


1999 






2155-304 


XBL 


0.097 


12.03 


10.88 


1979- 


1999 


0.66 


W92 


2200+420 


RBL 


-0.019 


4.63 


2.72 


1968- 


1999 


-0.05 


F03 


2254+074 


RBL 





8.24 


4.88 


1979- 


1999 


1.77 


W92 



Note to the Table: Col. 1 gives the name of the source, Col. 2 the classification, XBL is for X-ray selected BL 
Lacertae object, RBL for radio selected BL Lacertae object while mid-term for the object between XBL and RBL, 
Col. 3 for the spectral index(iS„ oc u~ a ), Col. 4 for the averaged observation polarization at 8GHz, Col. 5 for the 1 a 
deviation of the averaged polarization in Col. 4, Col. 6 for the time span of the observation data, Col. 7 for the 
core-dominance parameter, Col. 8 reference for the core-dominance parameter in Col. 7, respectively. 
F03: Fan & Zhang (2003); G93: Ghisellini et al. (1993); W92: Wills et al. (1992) 



